Spontaneous toroidic effects in Ba2CoGe207 
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The unusual magnetoelectric effects observed in the multiferroic phase arising below T/v=6.7K 
in Ba2CoGe2 07 (BCG) are related to the spontaneous toroidal moment existing in this compound. 
The transition to the multiferroic state, which involves spontaneous magnetization, polarization and 
toroidal moment gives rise to spontaneous toroidic effects. These effects produce specific contribu- 
tions to the spontaneous polarization and magnetization under applied magnetic or electric fields 
which provide indirect indications of the existence and role of the toroidal moment in multiferroic 
materials. The toroidic contribution to the electric polarization in BCG is shown to result from 
single-ion effects. 

PACS numbers: 77.80.-e, 61.50.Ah, 75.80.+q 



The resurgence of interest in multiferroic materials 
has prompted discussion of the relevance of the concept 
of magnetic toroidal moment for clarifying the macro- 
scopic and microscopic properties of these systems. 
The existence of a macroscopic moment asymmetric un- 
der both time reversal and space inversion long remained 
elusive^^^ until the observation of the independent co- 
existence of ferrotoroidic and antiferromagnetic domains 
in LiCoP04J^ This result provides a motivation for in- 
vestigating toroidic effects in the ferroelectric phases 
of magnetic multiferroic materials, in which the space- 
asymmetric electric polarization is induced by a time- 
asymmetric and space-asymmetric magnetic order. In 
absence of well defined physical properties showing direct 
experimental evidences of a toroidal moment in magnetic 
systems, one of the important issues is to find a material 
in which specific magnetoelectric effects^ would reflect 
indirectly the influence of the toroidic moment. Here 
we analyze theoretically the magnetoelectric effects dis- 
closed in the multiferroic phase of Ba2CoGe207 (BCG)^^ 
and show that due to its specific magnetic symmetry, 
which allows existence of spontaneous magnetization, po- 
larization and toroidal moment as well as linear and non 
linear magnetoelectric effects, the existence of a spon- 
taneous toroidal moment gives rise to spontaneous and 
field-induced toroidic effects. At the microscopic level the 
toroidic contribution to the electric polarization is shown 
to result from single-ion effects. 

The P42imV paramagnetic space group of BCG has 
at the centre of the tetragonal Brillouin-zone (/c=0) five 
irreducible representations (IRs) denoted ri — r^.^^ The 
one-dimensional IR's ri — induce non-polar magnetic 
symmetries P42im, P42im', Pi'2im and P4'2im', and 
are not associated with the transition to the ferroelec- 
tric phase observed in BCG below T/v =6.7K. The 2- 
dimensional IR rs spanned by the order-parameter com- 
ponents r]i = pcos{0) and r]2 = psin{6) is associated with 
the Landau expansion: 



Minimizing F yields 3 possible magnetically ordered 
phases below the paramagnetic phase: 



and 



1) Phase I corresponds to = ± 

cos {40) = 1 {0=n^). It has the orthorhombic mag- 
netic space groups P2;2i2' {r]i = p^7^2 = 0) or P2i2;2' 
(^1 = 0,772 = p^), which form energetically equivalent do- 
mains of the same equilibrium phase. 

2) Phase II has the magnetic symmetry Cm'a2' with 
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and cos{4.0) = -1 {0 = (2n + l)f ), 

or equivalently r]i = ±772 = p^- The phase involves a 
variety of spontaneous physical properties and domains 
represented in Fig. [H Denoting si and 52 the magnetic 
spins associated with the Co^+ ions located at (0,0,0) and 
(0.5,0.5,0) positions, the phase displays four weak ferro- 
magnetic/antiferromagnetic domains with a spontaneous 
unit cell magnetization M = si -\- S2 along the Tn'-plane 
and unit cell antiferromagnetic vector L = si — S2 along 
the a-plane. The polar symmetry of the phase also gives 
rise to two ferroelectric (iP^) - ferroelastic {^Cxy) do- 
mains, and four toroidic domains corresponding to the 
spontaneous toroidal moment T = z>(M x P), collinear 
to L, where z> is a third rank tensor. 
3) Phase III of symmetry P2' is stabilized for ^ rij 
and requires an eighth-degree expansion of P, involving 
8 weak-ferromagnetic/antiferromagnetic domains and 2 
ferroelectric- ferroelastic domains. 

The theoretical phase diagram containing the preced- 
ing phases is shown in Fig. [2l Phases II and III allow a 
spontaneous polarization along as observed experimen- 
tally in BCG^^^. However, only phase II can be reached 
directly from the paramagnetic phase whereas a transi- 
tion to phase HI goes across phases I or II, or displays 
a first-order character. The temperature dependence of 
the polarization varies continuously at Tat, and the con- 
figuration of the spin moments^ is consistent with the 
four domains of phase II. Therefore, phase II can unam- 
biguously be identified as the ferroelectric phase arising 
below Tat. In this phase F can be written in function of 




FIG. 1: (Color online) Respective orientations of the magnetization (M), antiferromagnetic vector (L), toroidal moment (T) 
and polarization (P) in the four multiferroic domains of BCG. Dark blue thin arrows represent the orientation of the spins 
si and S2 located in positions (0,0,0) and (0.5,0,5,0) in the tetragonal paramagnetic structure. The central inset shows the 
tetragonal and orthorhombic settings used in the text. 



M = (M^, My), L = {L^, Ly) and as: 

F = aiL^ + a2L^ + hM^ + b2M^ 

p2 

c{M^L^ - MyLy) + -j- + 5lL^LyP^ 



52M^MyP, + SsiM^Ly - MyL^)P, 



(2) 



where e^^ is the dielectric permittivity in the paramag- 
netic phase, and ai^bi^c and Si are phenomenological co- 
efficients. The equilibrium polarization below T/v reads: 

Pz = -4z i^iLxLy + S2M^My + SsiM^^Ly - MyL^)] (3) 

The two first terms into brackets express the re- 
spective contributions of the antiferromagnetic and 
weak-ferromagnetic order-parameters to the polarization, 
whereas the third term reflects their coupling contribu- 



tion. Since 



and 



M 



vary below T^v as ~ {T^ — T)^^ 

P| varies as (Tat — T), consistent with the linear depen- 
dence observed for P^{T)^^ corresponding to an improper 
ferroelectric critical behaviour. The dielectric permit- 
tivity varies as ezz{T) — e^^ ^1 + ^fSf^) ^ ^ 
in agreement with the reported upward discontinuity at 

The (^a-term in Eq. (|2]) reflects the invariance of the 



mixed vector product {M x L) • P under the symme- 
try operations of the paramagnetic phase. Analogously, 
the existence of the {T^My — TyMx)Pz invariant involv- 
ing the toroidal moment components {T^^Ty)^ expresses 
the invariance of the mixed vector products (T x M) • P 
and (T X P) • M which yield the following relationships 
between the spontaneous components: 



and 



P = jl{Tx M) 



M = X{T X P) 



(4) 
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where jl and A are third rank tensors. Applying electric 
or magnetic fields the existence of a spontaneous toroidal 
moment gives: 



P = P'' +eE + aH + a" {f x H) 



M = + XH + PE + a^{T x E) 



k^E 



d^"{E X H) 



(6) 
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where the third-rank tensors and cr^ precede ad- 
ditional polarization and magnetization components in- 
duced by the coupling of the total toroidal moment T 
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FIG. 2: (Color online) Theoretical phase diagram associated 
with the free-energy F given by Eg. ([2]). Solid and hatched 
curves are first and second-order transitions lines. Hatched- 
dotted curves are limits of stability lines. Ti and T2 are tri- 
critical points. N is a four-phase point. The arrow represents 
the thermodynamic path followed in BCG. 

to non-collinear H or E fields, and are the elec- 
trotoroidal and magnetotoroidal tensors. The (j^^-term 
represents the induced toroidal contribution under non- 
collinear electric and magnetic fields. One should em- 
phasize that the toroidal contributions to P and M in 
Eqs. (j6j) and ([7]), which are activated by a single exter- 
nal field non-collinear to T, differ from the " ferromagne- 
totoroidic" and " ferroelectrotoroidic" effects^^ implying 
the application of two distinct non-collinear fields. Al- 
though these contributions represent higher-order effects 
they can be differentiated from the non-linear magne- 
toelectric contributions. For example, the non-linear 
contribution to P is observable above and below the tran- 
sition while the T x H contribution is observable only 
below. More generally, the symmetries of the third-rank 
tensors and are different from the symmetries of 
the non-linear magnetic and electric susceptibilities, and 
the critical behaviour at constant fields of the correspond- 
ing tensor components is also different. 

Eq.(|6]-to-(j8])) provide an interpretation of the remark- 
able magnetoelectric effects reported in BCG.«^'^ Apply- 
ing Hz field gives rise to a polarization component Px 
increasing from to 120/i(7m~^ for < Hz < ST. In 
the orthorhombic setting, which is turned by 45° with 
respect to the tetragonal axes (Fig. [T]), one gets from 
Eq. dHj): 



P.iHz) = [ai3 + T^ixss + af23)] H 



(9) 



Consistent with the linear increase observed for Px{Hz) 
with increasing field^ and with its sign reversal on re- 
versing Hz.^ The sharp increase of PxiT) below Tn at 
constant field, ^ denotes a substantial toroidal contribu- 



tion (J 



H rps 



(Tat — T) 2 , with respect to the linear 
magnetoelectric contribution 0^13 ^ (T/v — T)^. The ob- 



served increase of Pz{Hz) from — 11/iCm ^ at iJ^ = 0, 
to +80/iCm~^ in 8T is given by: 



(10) 



The shift of Pz{Hz) to higher temperature under ap- 
plied field^ is due to the renormalization of the coeffi- 
cient ai ^ {T — Tn) in Eq. (|2]), which increases Tn 
by Tn{Hz) — Tn{0) ^ X33^2- order to account for 
the even dependence of Hz observed for Pz{Hz)^ one has 
to consider a higher order contribution, e.g. :^ H^^ to 

Other magnetoelectric effects have been reported^ un- 
der application of Hxy and Hxy fields. Pz increases by 
increasing Hxy and decreases when increasing Hxy.^ Pro- 
jecting Eq. ® along one gets: 



Pz{Hxy) — —APz{Hxy) = 

-(^31 + (T32lTy + ^312Tx)Hi 



(11) 



Turning the Hxy field by 90° transforms a ferroelectric 
domain into another, changing the sign of P^.^ As for 
Pz{Hz), a shifting of the transition temperature is ob- 
served under Hxy field^ Pz{T) decreasing smoothly down 
to TAr=12K for Hxy=hT , with a {Tn - T)i critical de- 
pendence of A{T^,T^) = (J^2iTy + ^3i2Tx' These effects 
occur at low magnetic fields. At higher fields Px{Hxy) 
decreases and changes sign.^ This behaviour, assumed 
to correspond to a spin-structural change,^ requires in- 
cluding the higher-order invariant (T^ x H) • (T* • H) ^ 

^^^^^H^y in Eq. For K < 0, Pz{Hxy) decreases 
above the threshold field H^y = — ^^^^"^ taking negative 
values for Hxy > 2i^^^. 

To gain insight into the nature of the magnetic inter- 
actions governing the magnetoelectric and toroidic be- 
haviours of BCG, let us express the order-parameter 
components in function of the magnetic spins si and 
52- Writing si = sfa + s^b -h s^c {i = 1,2), where 
a, 6, c are the tetragonal lattice vectors, the representa- 
tion r transforming the 5^'^'* components decomposes 
into r = Ti + r2 + 2r5, i.e. two order-parameter copies, 
denoted (7^1,7^2) and ctre involved in the transi- 

tion mechanism. Standard projector techniques'^ give: 



Ci 



S2,V2 = - 
'1 ~ ^25 S2 
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It shows that the two order-parameter copies coincide 
with the ferromagnetic and antiferromagnetic vectors. 
On the other hand, projections of T on ri and r2 lead 
to 5^ — ^2 = and s^; + 52 = 0, i.e. s^ = S2 = 0, confirm- 
ing the in-plane spin ordering in BCG. The equilibrium 
values of (7^1,7^2) and {C1X2) in phase II yield the spin 
configurations for the four magnetic domains represented 
in Fig 1, namely: two weak ferromagnetic domains for 
s^-\-S2 = ±(5^+52), and two antiferromagnetic domains 
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for s\ — S2 = ^{si — S2) . The spontaneous polarization 
P| at zero field reads: 

= s'lmm + 52C1C2 + ^3(^(2 + mCi) (i3) 

analogue to Eq. (3). Using Eq. (12) yields: 

+5',isf - sf - sf + sf ) (14) 

Eq. ([H]) holds for a pair of antiferromagnetic domains 
(e.g. r]i = 772 and = Q) whereas — P| coincides 
with the other pair {r]f = — r^f^Ci = "CI)- The 
62 and terms represent the respective contributions 
of the spontaneous ferromagnetic, antiferromagnetic and 
toroidal contributions to the spontaneous polarization 
arising in the multiferroic state. The ^3 term, which is the 
microscopic analogue of the spontaneous toroidic effect 
given by Eq. (|4]), reflects single-ion effects^ while the two 
other terms contain invariants sfs^{i = 1^2] u^v = a, 6), 
also corresponding to single- ion effects, and sfs^{i ^ j) 
invariants expressing the symmetric part of the exchange 



coupling interaction between the two Co spins. These re- 
sults support the interpretation^ that the spin-dependent 
p—d hybridization between the transition- metal (Co) and 
ligand (O) contributes to the ferroelectricity in BCG via 
the spin-orbit interaction, as well as the proposed mecha- 
nism of lattice relaxation induced by exchange striction^^ 
Note that the Dzialoshinskii-Moriya (DM) interaction 
does not contribute directly to the polarization but is re- 
sponsible of the canting inducing the weak-ferromagnetic 
moments^ which stabilizes the toroidal moment giving 
rise to the ^3 term in Eq. (p!4|) . 

In conclusion, our theoretical analysis shows that in 
multiferroic compound exhibiting spontaneous magneti- 
zation, polarization and toroidal moment, field-induced 
toroidic effects occur, consisting of additional toroidal 
contributions to the polarization and magnetization. 
These toroidal effects allow a comprehensive descrip- 
tion of the specific magnetoelectric properties observed 
in BCG, which have been shown to relate to additional 
toroidal contributions to the polarization, corresponding 
at the microscopic level to single site magnetic interac- 
tions. They should allow clarifying the intrinsic role of 
the toroidal moment in magnetic multiferroics. 
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